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ABSTRACT 

We have detected a strong deflection of radio waves from the pulsar PSR B0834+06 in scintillation 
observations. Interference between the undeflected pulsar image and deflected subimages allows single 
dish interferometry of the interstellar medium with sub-milliarcsecond resolution. We infer the pres- 
ence of scattering structure(s) similar to those that are thought to cause Extreme Scattering Events 
in quasar flux monitoring programs: size ^ 0.2 AU (an angular size of 0.1 mas) with an electron 
overdensity of > 10^ compared to the warm ionized medium. The deflectors are nearly stationary in 
a scattering screen that is thin (< 5% of the pulsar-observer distance in extent), is located 70% of the 
way from the Earth to the pulsar, and has been seen consistently in observations dating back 20 years. 
The pulsar scans the scattering screen at a velocity of 110 km s~^ with a detection radius of 15 mas. 
Pulsar observations such as these — particularly with a new generation of low-frequency radio tele- 
scopes with large collecting areas — hold promise for improving constraints on the poorly understood 
physical characteristics and space density of the deflecting structures. Such observations may also 
prove useful in correcting deviations the deflectors produce in high-precision timing of millisecond 
pulsars. 

Subject headings: ISM: structure — pulsars: general — scattering — techniques: spectroscopic — 
pulsars: individual (PSR B0834-I-06) 



1. INTRODUCTION 

The discovery o f Extreme Scattering Events (ESEs; 
iFiedler et alJll987^ in the late 1980s provided evidence 
for ultracompact (< 1 AU), ionized objects in the inter- 
stellar medium. Despite more than 15 years of obser- 
vational and interpretive work, the physical conditions, 
origin, and lifetime of these compact refractors are not 
well establish ed (Romani et al. 1987"; "Fiedler ct al. 1994'; 
IRickett et a].iri997t Wa,1ker Wa,rdle. ,1 998^. Further- 
more, their space density and possible contribution to 
the mass distribution in the Galaxy are poorly known. 
The basic difficulty has been detectin g events in ei ther 
quasar or pulsar monitoring programs l)Walkeill2001[) . In 
order to substantially affect the flux density of a radio 
source, the refractor must be nearly aligned with the 
source. The refractors are typically located at distances 
of kiloparsecs, so a milliarcsecond alignment is required. 
Thus, only about a dozen ESEs have been identified, and 
few of these are easy to model. 

The requirement of stringent alignment is relaxed if 
the intrinsic source size is small enough to exhibit in- 
terstellar scintillation, one of the observational conse- 
quences of coherent multipath scattering from density 
inhomogeneities in the medium. Pulsars are the domi- 
nant source in this category. Interference between var- 
ious parts of the pulsar image gives rise to faint but 
observable periodic fringing patterns in pulsar dynamic 
spectra (radio flux density as a function of time and ob- 
serving frequency) . These features have a clear represen- 
tation in the secondary spectrum (the squared modulus 
of the Fourier transform of the dynamic spectrum) and 
yield information about the location of the scattering ma- 
terial and the structure of the image (Roberts & Ablest 
[1982: Hcwish et al. 1985: Cordes fc YiQlmcz m 19S 
IWolszczan fc CordedligSTt IRickett et al.lll997f) . 
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ticular, periodic fringing in the dvn amic fring ing gives 
rise to parabolic scintillation a rcs ( Stine^m^^ al. 2002 
IHill et all200^ICordes et al.l2004i:iWalker et alJ20 04) in 

the secondary spectrum, which imply that the dominant 
scattering material is localized in a thin screen along the 
line of sight. 

In this Letter, we present evidence for multiple imag- 
ing of the pulsar PSR B0834-t-06 with greater detail than 
a similar event reported by Rickctt ct al. (1997). We de- 
tected secondary spectrum features caused by at least 
five separate ray paths throughout a 26 day observing 
run. The angular separation between the undefiected 
pulsar image and each of the four deflected subimages 
grew linearly at a rate consistent with the velocity of the 
pulsar. This, combined with multi-frequency data, im- 
ply that the multiple imaging is caused by one or several 
deflecting structures similar to those that cause ESEs 
and that the structures are stationary in the scattering 
screen. S uch objects, if they trace out underlying neutral 
material ijWalker fc War dle' 1998), may contain a signif- 
icant fraction of the mass of the Galaxy. 

2. THIN SCREEN MODEL 

We use a model of scintillatio n arcs developed in 
IStinebring et~an IHifl et all ll200l . ICordes et alJ 

Ij2004ll . a,Tid IWalker et al.l l)2004Jl and summarized here. 
The essential feature of the model is that a scintillation 
arc is caused by scattering in a thin screen along the 
line of sight to the pulsar. Interference between a bright 
core of the pulsar image and the scatter-broadened halo 
gives rise to the basic arc, and the details of the power 
distribution along the arc provide information about the 
power distribution in the pulsar image. Therefore, our 
secondary spectra of PSR B0834+06 at 327 MHz provide 
a partial image of the pulsar with an angular resolution 
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of ~ 0.3 mas and a field of view of ^ 40 mas (see below). 

2.1. Image features and interference patterns 

We define a coordinate system in the plane of the sky 
with the pulsar at the origin, the O^-sods pointing along 
the effective velo city vector, and 6,, -L 6^.. The pulsar 
effective velocity l|Cordes fc Ricketti (l99S1. 
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Veff = (f - s)V^ 



sVobs - Vs 



(1) 



is the apparent velocity of the image of the pulsar 
through the thin scattering screen, where s is the frac- 
tional position of the screen from the pulsar (s = 0) to 
the observer (s = 1), Vp is the transverse velocity of the 
pulsar, Vobs is the transverse velocity of the observer, 
and Vs is the transverse velocity of the screen. Interfer- 
ence between two arbitrary points in the image plane, 
6i and 02, causes fringing in the dynamic spectrum 
with conjugate time and conjugate frequency of ft = 
-{62 - 9,) ■ Feff/(sA)^ and = D{1 - s){el - el)/{2sc), 
where A is the observing wavelength. The fi, coordinate 
is a measure of differential time delay between pairs of 
rays, and ft represents the temporal fringe frequency of 
the interference between the two rays or, alternatively, 
the differential Doppler shift between them. Interfer- 
ence between the ray at the origin and points along the 
^X-axis produces a parabolic scintillation arc defined by 
fv — vft I where the arc curvature parameter is 



V ■ 



DX^s{l - s) 
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Interference between the origin and points with non-zero 
6y places power inside the parabola because (x 9'^+9y. 
Interference between a bright spot in the periphery of the 
image and the the rest of the image produces an inverted 
parabola or arclet with the same I77I and a vertex with 
an ft coordinate that is related to the 0a;-coordinate^ of 
the bright spot by 



(3) 



Such features are accentuated when the image is elon- 
gated along the velocity vector or the bright spot lies 
near the ^^j-axis, or both. 

2.2. Screen Location 

For pulsars with measured proper motions and mea- 
sured or estimated distances, we can determine the lo- 
cation of the dominant scattering material from the 
curvature of the main scintillation arc using equa- 
tion O . PSR B0834+06 has an estimated distance 
ijCordes fc Lazia2003rl and a mod eratelv well determined 
proper motion l)Lvne et"allll982() : D = 0.64 ± 0.08 kpc 
and /i = 51 ± 3 mas yr""'^. We have measurements of 
scintillation arcs in PSR B0834+06 dating back to 1981 
(Stinebring et al. 2005, in preparation). Those data 
and the observations presented here are consistent with 
a value of 77 = 0.47 ± 0.03 s^ at 327 MHz. This resuhs 

^ The minus sign in equation El and in the defintion of ft is 
present because negative ft values correspond to deflected rays in 
front of the moving pulsar; points behind the moving pulsar (past 
closest encounter to the deflecting structure) produce a positive ft 
feature. 
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Fig. 1. — The dynamic spectrum of PSR B0834-I-06 observed on 
2003 Dec 31. The flux density as a function of frequency and time 
is shown using a grayscale that is linear in power, with dark regions 
indicating high power. The crisscross pattern is due to radio waves 
reaching the observer from a variety of angles (~ 10 mas away from 
the pulsar position), as detailed in the text. 



in a value of s = 0.29 ± 0.04 or a screen that is located 
a distance 0.46 ± 0.08 kpc from the observer. These pa- 
rameters result in a conversion between /( and angle on 
the sky of 50 mHz = 24 mas. Hence, the angular reso- 
lution of the secondary spectrum is (2 x 24 mas / 180 
pixels) — 0.27 mas per pixel; the useful field of view is 
about 40 mas. 

3. OBSERVATIONS 

To search for evidence of compact refractors using pul- 
sar scintillation, we observed the pulsar PSR B0834-I-06 
at the Arecibo Ob servatory i n 2004 January using the 
same technique as IHill et aT ? f2003'). An example dy- 
namic spectrum is shown in Figure ^ We obtained dy- 
namic spectra during 30-60 minute integrations using 
the 327 MHz receiver with spectral resolutions of ap- 
proximately 1.5 kHz on 11 days over a 26 day period. 
Typically, we simultaneously took data with center fre- 
quencies of 321 and 334 MHz using multiple Wideband 
Arecibo Pulsar Processor spectrometers. The secondary 
spectrum in Figure|21exhibits a scintillation arc as well as 
numerous arclets. The four isolated arclets, labeled a-d, 
shifted upward and to the right along the main parabola 
during the month. This had never been seen before and 
is the central observational result of this paper. 

4. MOTION OF ARCLETS 

We developed an algorithm to locate the position of 
these arclets in secondary spectra by summing along 
parabolic regions with vertices incremented along the 
main parabola. In Figure |31 we plot the fringe frequency 
of arclet a as a function of day number for two closely 
spaced frequencies. Two models of the refracting struc- 
tures predict different frequency behavior: if the rays are 
refracted by smoothly varying electron density variations 
, the angular position of the rays would sc ale as is typica l 
for the interstellar medium: 6* oc (e. g. lScheuerlll968(l : 
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Fig. 2. — The secondary (or delay-fringe-frequency) spectrum 
corresponding to the dynamic spectrum in Figure^ The grayscale 
is logarithmic in power and represents a total range of alsout 10^ 
from the maximum power to the noise level. A primary arc with 
curvature r; = 0.47 s"^, fitted to the left side of the main parabola, is 
shown with a dashed line. Numerous inverted parabolas, or arclets, 
are present with vertices near the primary arc. The four isolated 
arclets (labeled a-d) each correspond to a distinct enhancement 
in the pulsar image. The horizontal axis can be converted into 
a separation angle projected along the direction of pulsar motion 
using 2 mHz fs; 1 mas. 



then, ft oc d\~^ oc A. Alternatively, if the rays giving 
rise to the arclets pass through a discrete an d dense re- 
fracting structure (lens model) ijClegg et al.l ll998). they 
would arrive from nearly the same angle at different 
frequencies; they still deflect as within the structure, 
but only an unmeasurably slight deviation with wave- 
length occurs because of the small size and high density 
of the refracting structure. Thus, ft oc dX^^ oc A^^. 
We performed a least-squares linear fit to the lower fre- 
quency data. We then calculated lines with the same 
slope but vertical offsets scaling as in the lens model and 
in the smooth variation model. The data in Figure |31 
clearly support the lens model, which we assume in the 
following discussion. 

Using measured values of ft for the vertices of the ar- 
clets, we calculated their angular position as a function 
of time, as shown in Figure 01 Uncertainties in the an- 
gular positions were typically 0.2 mas as derived from 
the scatter about the best fit lines. Arclet a moved at 
a rate of 47 ± 2 mas yr~^ and arclets b-d each moved 
at 51 ± 2 mas yr~^ along the pulsar effective velocity 
vector. Taking the measured proper motion of the pul- 
sar and accounting for the Earth's motion a t this epoch 
l|van de KamiJl967tlCordes fc Rickettll998|) yields an ef- 
fective pulsar angular velocity of /Zeff = 50±3 mas yr^^. 
Hence the angular motion of the pulsar, determined 
from astrometry, and the rate of increase of the pulsar- 
deflector angle, determined from the data reported here, 
are consistent within their combined uncertainties. We 
use this comparison to estimate the motion of the deflec- 
tor relative to the screen. We find Wgcr = (1 — •s)-D(/Xarc — 
fJ-cs) ^ 9 km s~ , implying that the physical objects giv- 
ing rise to the arclets are essentially stationary in the 
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Fig. 3. — The fringe frequency (ft) coordinate of the vertex of 
arclet a versus day number (day = 2003 Dec 31). Data taken at 
321 MHz are plotted with filled circles; 334 MHz data are plotted 
with open circles, and the uncertainty in position is about the 
diameter of the symbols. The solid line represents a two parameter 
fit to the 321 MHz data. The dashed line has the same slope, 
but with the offset expected if /t oc A~^, or 6 is constant with 
observing wavelength; the dotted line has the offset expected if 
/i oc A, or scattering angle 9 oc These data indicate that a 
lens-like structure or structures in the scattering screen causes the 
arclets. 
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Fig. 4. — The angular position of four arclets over 26 days of 
observations at two observing frequencies (day = 2003 Dec 31). 
As in Figure|21 filled circles denote 321 MHz data and open circles 
indicate 334 MHz data. The uniform, consistent motion of the 
arclets is primarily due to the motion of the pulsar, indicating that 
arclets are caused by scattering objects essentially stationary in 
the screen. 



scattering screen. 



5. DISCUSSION 



The refracting angle of a plasma lens of size a is 
9r ~ aa/D, where a = X^r^NQD /no? characterizes the 
strength of the refractive effects, Nq = J n^ds is the 
maximum electron column density through the lens, and 
re is the classical electron radius (jClegg et al.lll998() . If 
we define an average electron density fie = Nq/q, then 
fie — (5.4 cm~^)9r/X'^, where the angle is measured in 
mas and the wavelength is expressed in meters. For 
9r « 15 mas at 327 MHz, we find fie ~ 100 cm~'^, which 
is more than lO'^ times greater than t he average elec- 
tron density in the interstellar medium l|Cordes fc Laziol 
I2OOI . 

The thinness of the arclets is remarkable as well and re- 
inforces the interpretation of refraction by lens-like struc- 
tures embedded in a thin scattering screen. The arclets 
in Figure |2] have a fractional thickness of A/^//^ « 5%. 
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Since oc 9'^, this implies a - A6'(l - s)D w 0.1 AU if 
each arclet is produced by a separate refractor. In this 
case the mass in each refractor is only M ~ a^fiemH ~ 
10~i^ Mq; here mn is the mass of a hydrogen atom and 
we are assuming a fully ionized object. 

Instead of four separate deflecting structures, the four 
arclets may be caused by two or even just one struc- 
ture. The basic imaging features of a plasma lens be- 
come apparent by analyzin g a sirnple m odel with a Gaus- 
sian column density ijClegg et al.P l99^ . This model pre- 
dicts caustic surfaces where rays intersect and a defo- 
cusing region directly in front of the lens. In the region 
between the inner caustic and the outer caustic an ob- 
server would see a bright undeflected image of the source 
and two dim, deflected subimages, which would give rise 
to two arclets as the subimages interfere with the cen- 
tral image. The angular separation between a pair of 
arclets remains constant, as in these observations, for 
the arclet pairs ab and cd, if the refractive strength pa- 
rameter a ;2> 1; in this case, A9 « 2a/{D — sD) and 
fie = na{l ~ s)A9/{2X^re) « (1.9 cm-3)aA0„ias/ALtor- 
We can estimate a from the observations because the 
brightness of the deflected image relative to the central 
image is « {1 + a)^^. The peak power of the arclets 
in Figure 13 is about 10"'^ of the power at the origin in 
the secondary spectrum, so a > 10'^. Using this model, 
our observations yield estimates of a « 0.25 AU, fie ~ 
2000 cm-3, and M « 4 X IO-^'^Mq. It should be noted 
that neither of these models is in pressure balance in the 
warm ionized phase of the interstellar medium; hence, 
they would be expected to dissipate quickly unles s con- 
strained through some other mechanism ( Roman i et aD 
119871 mickett et al.l[T997t [ Walker fc Wardl&.199^1 . 

IWalker &: Wardld l|1998|) proposed a model for ESEs 
in which compact (~ 1 AU) , long-lived spheroids of neu- 
tral hydrogen are surrounded by a photoionized shell. 
This results in a double-peaked electron column density 
profile; each peak may give rise to two enhancements 
in the image, resulting in the observed four-arclet pat- 
tern. In this case, a « A9ad{^ — s)D/2 « 1 AU, where 
A9ad = 3.3 mas is the separation between arclets a and 
d. The mass encompassed by the refracting structure is 
dramatically higher in this model (M ~ 1O~^M0) be- 
cause of the unionized, cold, dense core. A Galactic halo 
population of these clouds, sufficient to account for the 
poorly constrained ESE rate, would be a major contrib- 
utor to the total mass of the Galaxy. 

Because interference effects are present even when the 
deflector does not intercept the line of sight, pulsar scin- 
tillation observations will improve constraints on the 
space density of the deflecting structures. The angular 
offset of the deflector(s) from the 0a;-axis is an important 



parameter in estimating their space density. Analysis of 
these data yields 9y offsets in the range 2-6 mas. Our 
sensitivity to this level of offset, along with the ability to 
detect deflectors within 10-15 mas on either side of the 
pulsar, greatly increases the detection probability com- 
pared to flux-only observations. 

6. CONCLUSIONS 

Scintillation observations provide greater detail on the 
physical properties of these refractors than quasar flux 
monitoring for two principal reasons. First, the inter- 
ference between the deflected and undeflected images of 
the pulsar allows a single dish telescope to act as an 
interferometer with excellent angular resolution. Sec- 
ond, pulsars scan the scattering screen quickly because of 
their high transverse velocity. The data presented here 
demonstrate that dense, compact refractors similar to 
those that give rise to ESEs can be detected in pulsar 
scintillation observations. The refracting structures are 
of interest in their own righ t, but they also affec t high- 
precision timing of pulsars (Cog nard et alJlT99^ . Since 
this is one of the most promising methods of detecting a 
background of primordial gravitational waves, the abil- 
ity to detect - and potentially correct for - the effects of 
refracting structures is of practical importance as well. 

We currently do not have enough data containing 
arclets to decide between models or to place realis- 
tic limits on the space density of refracting struc- 
tures. Arclets are not rare, however, particularly at 
lower frequency (< 400 MHz) where the effects of 
scattering and refraction are more pronounced. We 
have seen arclets, on occasion, in data from pulsars 
(PSR) B0355+54, B0823-h26, B0834-h0 6, B1133-H6. 
B1642-03, B1737-hl3, and B1919-h21. iRickett et ail 
( 1997) identified what is probably a poorly- r esolve d ar- 
clet in B0834+06, and Gor des fc Wolszcza^ lll986'^ and 
Wolszczan & Cordcs ( 1987) report similar features in ob- 
servations of B0919-h06 and B1237-h25. We are analyz- 
ing all of these data further in order to obtain a realistic 
estimate of the space density of deflectors. It is clear, 
however, that single-dish scintillation observations with 
adequate sensitivity at low frequency can be used to track 
deflectors with milliarcsecond resolution and to estimate 
their physical properties. 
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